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The reconstruction of quantum states from a sufficient set of experimental data can be achieved 
with arbitrarily weak measurement interactions. Since such weak measurements have negligible 
back-action, the quantum state reconstruction is also valid for the post-selected sub-ensembles usu- 
ally considered in weak measurement paradoxes. It is shown that post-selection can then be iden- 
tified with a statistical decomposition of the initial density matrix into transient density matrices 
conditioned by the anticipated measurement outcomes. This result indicates that it is possible to 
ascribe the properties determined by the final measurement outcome to each individual quantum 
system before the measurement has taken place. The "collapse" of the pure state wavefunction in a 
measurement can then be understood in terms of the classical "collapse" of a probability distribution 
as new information becomes available. 
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I. INTRODUCTION 

As our ability to control individual quantum systems 
increases, the seemingly paradoxical aspects of quantum 
measurement theory take on a more practical relevance. 
One recent example is the renewed interest in the quan- 
tum statistics of post-selected weak measurements, which 
seem to suggest the presence of negative probabilities as 
the source of quantum paradoxes El 0, EL IE H> S Hi • 
These developments in the field of weak measurement 
may be especially significant in the context of new exper- 
imental possibilities pioneered by quantum information 
related research @, [13, UH- However, there seems to be 
a certain mismatch between the conventional approach 
to weak measurements, which is based on the notion of 
measurement interaction dynamics mediated by opera- 
tor observables [H, M M M M , 

and the more general 
approach to measurements based on the measurement 
operators widely used in quantum information [l7| . In 
particular, the conventional analysis seems to overem- 
phasize the exotic and surprising aspects of the weak 
measurements, while the operator based approach shows 
more clearly how weak measurements fit into the gen- 
eral framework of quantum physics [l8| . As experimental 
weak measurements become more and more established, 
it may therefore be time to shift the focus away from the 
oddities of specific cases, towards a complete and consis- 
tent formulation of post-selection effects in terms of their 
experimentally observable properties. 

From the experimental side, quantum states and pro- 
cesses can be characterized by measuring their complete 
statistical properties, a procedure known as quantum to- 
mography |l9l . [20| . A significant merit of quantum to- 
mography is that it establishes an operational approach 
to quantum states, that is, it defines the quantum state 
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in terms of the experimentally accessible data. By apply- 
ing quantum tomography to generalized weak measure- 
ments, it is possible to extend this operational definition 
to post-selected sub-ensembles of a quantum state. In 
the following, the general theory of quantum tomography 
with measurements of variable strength is formulated. 
It is shown that, in the limit of weak measurements, 
post-selection partitions the initial density matrix into 
sub-ensembles described by non-positive transient den- 
sity matrices. This result divides the "collapse" of the 
wavefunction into two distinct parts, one associated with 
the selection of the appropriate sub-ensemble, and the 
other related to the actual back-action of the measure- 
ment dynamics. Although the measurement back-action 
is necessary to cover up the negative eigenvalues of the 
transient states, weak measurement tomography suggests 
that the sub-ensemble partition has a physical meaning 
even before the measurement interaction takes place. 



II. QUANTUM STATE TOMOGRAPHY WITH 
WEAK MEASUREMENTS 

The starting point for the following derivation of a 
complete and consistent theory of weak measurement to- 
mography is the representation of general quantum mea- 
surements by a set of operators {M m } acting only on the 
Hilbert space of the system pjj ■ These operators summa- 
rize the relevant effects associated with a measurement 
outcome m, separating the essential properties of a quan- 
tum measurement from the technical problem of its im- 
plementation by a specific combination of system-meter 
interaction, meter preparation, and meter read-out. In 
the case of Hilbert space vectors representing pure states, 
the application of a measurement operator M m to a state 
vector changes both the length and the direction of that 
vector. The new direction of the state vector then repre- 
sents the output state after the measurement, while the 
squared length represents the probability p(m) of obtain- 
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ing the measurement outcome to. If the quantum state 
is expressed in terms of the density matrix p il the prob- 
ability p{rn) is given by a product trace, 

p{m) = Tr {M m piMl } = Tr {ik^M™ Pl } . (1) 

The set of squared operators {M^M m } is the positive 
operator-valued measure (POVM) of the measurement 
probabilities. Since all probabilities must add up to 1, 
the POVM fulfills the completeness relation 

Y,MlM m = l. (2) 

m 

The POVM formalism describes the general relation be- 
tween experimental data and the quantum state. Specif- 
ically, eq.([T]) shows that the measurement probabilities 
p(m) arc linear combinations of the density matrix ele- 
ments. If the set of relations given by eq.fl]) is invert- 
ible, the complete density matrix can be reconstructed 
from the available set of measurement probabilities. In 
a d-dimensional Hilbert space, quantum tomography can 
thus be performed using any combination of d 2 linearly 
independent measurement operators M^Mm. 

Whether a POVM is invertible and therefore suitable 
for quantum tomography does not depend on the preci- 
sion of the measurement. It is therefore possible to re- 
construct the quantum state from arbitrarily weak mea- 
surements. To illustrate this point, it is useful to for- 
mulate the inversion procedure for POVMs with variable 
strength e. The strength or weakness of a measurement 
can be quantified directly by the closeness of the mea- 
surement operators to multiples of the identity operator 
1. A convenient way of representing a variably measure- 
ment strength e in the formalism is 

MlM m = w m {\ + eS m ) 
with y^w m =—^— and w m S m = — j— . (3) 

m m 

The measurement probabilities in eq. |T]) can then be ex- 
pressed in terms of the expectation values of a set of self- 
adjoined operators S m that is independent of the mea- 
surement strength e. For quantum state reconstruction, a 
set of d 2 linearly independent operators {S m } defines an 
operator expansion of the density matrix in terms of the 
set of reciprocal operators {A m } with Tr{S n A m } — S n>m . 
The density matrix is then given by 

Pi = 22Ti{S m pi}A m , (4) 

m 

where the coefficients of the expansion are related to the 
measurement probabilities p(m) by 

Tr{S m Pi} = — ^— (p(m) - w m ) . (5) 

ew m 



Eqs. (|4|) and (O express the density matrix in terms of 
the experimentally obtained measurement probabilities 
pirn). Thus, quantum tomography can provide a defini- 
tion of the density matrix that is based only on empirical 
properties of the system. 



III. QUANTUM STATES OF POST-SELECTED 
ENSEMBLES 

Conventional quantum tomography is a one-way read- 
out process in which the quantum state is discarded after 
the measurement. However, the measurement operators 
M m also provide a description of the quantum state after 
the measurement. It is therefore possible to consider the 
effects of a subsequent measurement with outcomes /, 
described by another POVM {II/}. The joint probabil- 
ity of obtaining first m and then / in the measurements 
is given by 

p(mj\i) =Tr{ft / M m piMl]. (6) 

Since the measurement operators and M m are not 
directly multiplied, the probability p(f\i) found by sum- 
ming over all to is different from the product trace of pi 
and Hf. This change in the probability of obtaining / 
from the initial state pi is caused by the measurement 
back-action associated with the measurement of m. It is 
therefore impossible to know whether the final result / 
was caused by the initial state of the system or by in- 
teraction effects related to the measurement outcome to. 
However, the theory of weak measurements shows how 
this problem can be circumvented: for very small mea- 
surement strengths e, the effects of the quantum state 
on the measurement probabilities is linear in e while the 
measurement back-action is quadratic in e. It is there- 
fore possible to realize quantum state tomography with 
negligible back-action. 

In the limit of weak measurements (e <C 1), the mea- 
surement operators M m are approximately given by the 
linearized square root of the POVM, 

M m « (l + J&n) • (7) 

In the joint probability p(m,f\i) given by eq.©, the 
terms linear in e are obtained by applying S m either 
from the right or from the left. If quadratic terms are 
neglected, this is equivalent to applying the POVM from 
the right or from the left. Therefore, the joint probabil- 
ity p(m, f\i) can be approximately expressed by the weak 
measurement POVM, 

p(m, f\i) w Tr jfty ^{p t MlM m + M^M^)} . (8) 

Since the POVM fulfills the completeness relation given 
by eq. ©, the final measurement probabilities p{f\i) are 
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not changed by the measurement of m and 

p(/|i) = X>( m '/K) = Mn/Pi}. (9) 

m 

Hence the measurement back-action is negligible and the 
measurement results m merely identify the state condi- 
tioned by both initial preparation and final measurement. 
The conditional probability can then be written in terms 
of a transient density matrix Ri t , such that 

P(m\i, /) = = Tr {MlM m Ri f } . (10) 

According to eqs.((8]) and ©, this transient density ma- 
trix is given by 

Since quantum tomography provides an unambiguous 
definition of the density matrix in terms of the mea- 
surement statistics of m for any well-defined ensemble of 
quantum systems, the consistency of quantum measure- 
ment theory requires that the self-adjoint operator Rif is 
the proper statistical representation of the sub-ensemble 
defined by the conditions i and /. 

IV. ANTICIPATORY DECOMPOSITION 

The operator Ri f correctly predicts the outcomes of all 
(real) weak values that can be obtained between i and /. 
In this sense, it is similar to the two-state vector formal- 
ism of weak measurements [2~l| and its extension to mixed 
states [22| . However, the derivation from tomography en- 
sures that Rif is a self- adjoined operator with real eigen- 
values and a trace of one. It thus corresponds to a con- 
ventional density matrix, except for the possibility of neg- 
ative eigenvalues. Moreover, the definition of Rif applies 
equally well to projective measurements of pure states 
and to noisy measurements of mixed states. The present 
analysis therefore bridges the gap between the classical 
limit and the extreme quantum limit, revealing similar- 
ities of quantum and classical statistics that tend to be 
obscured by the state vector formalism. Specifically, the 
transient density matrix Rif is the quantum mechanical 
representation of conditional probabilities p(m\i, /), just 
as the density matrix pi is the quantum mechanical rep- 
resentation of the probabilities p(m\i). The relation be- 
tween the total density matrix pi and the set of transient 
density matrices conditioned by the final measurement 
outcomes / is therefore given by the weighted sum over 
all possible outcomes /, 

Pi = $>(/|i) Rif- (12) 
/ 

This decomposition of the density matrix explains why 
the weak measurement statistics can be measured before 



the final measurement / has taken place: the measure- 
ment of / simply identifies a sub-ensemble Rif that is 
already included in the total ensemble pi. It is therefore 
possible to decompose pi into Rif in anticipation of the 
measurement of /. Such an anticipatory decomposition 
indicates that the measurement outcome / can already 
be ascribed to the quantum systems before the measure- 
ment has taken place. 



V. WAVEFUNCTION COLLAPSE WITHOUT 
MEASUREMENT BACK-ACTION 

In classical physics, the statistical rules for anticipatory 
decompositions are straightforward, since the state p is 
just a conventional probability distribution over all mi- 
croscopic configurations of the system. In the quantum 
case, it is usually assumed that the "collapse" of a pure 
state caused by a projective measurement is fundamen- 
tally different from such a probability update, because it 
includes the effects of decoherence. However, weak mea- 
surement tomography shows that the "collapse" can be 
divided into two parts, one associated with a classical 
sub-ensemble selection (that is, a Bayesian probability 
update), and the second one associated with the actual 
physical interaction that results in the decoherence. It 
is therefore possible to identify the set of quantum sys- 
tems that produce a specific measurement result / with a 
uniquely defined sub-ensemble of the total density matrix 
Pi- 

To understand the significance of this result, it may be 
useful to reconsider some of the paradoxes of quantum 
mechanics in the light of these experimentally accessi- 
ble facts. For instance, it is now possible to resolve the 
problem of double slit interference by assigning both an 
interference pattern and a slit to each individual par- 
ticle. Specifically, the initial pure state superposition 
| i) = (| 1)+ | 2))/y2 of a particle passing through slit 1 
and a particle passing through slit 2 can be decomposed 
into 

Ril = |l>(l|+i(|l)(2| + |2)(l|) 

Ri2 = |2)(2|+i(|l)(2| + |2)(l|) (13) 

in anticipation of a which-path measurement. Weak mea- 
surements performed between the preparation of the su- 
perposition | i) and the final which-path measurements 
show that the path information co-exists with the inter- 
ference pattern of the superposition. Therefore, the dou- 
ble slit interference pattern can be obtained from weak 
measurements even if the particle has passed through 
only one of the slits. Eq. (|13p thus shows that coher- 
ence between two alternatives does not imply that both 
alternatives are simultaneously realized. 

Essentially, weak measurement tomography is an anal- 
ysis of quantum statistics that reveals additional details 
about where quantum information is located before it is 
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measured. In particular, it fills a gap left by the measure- 
ment postulate by showing that the measurement result 
/ corresponds to properties of the system before the mea- 
surement, and is not just randomly generated by a role of 
the dice at the time of measurement. Thus, we can solve 
the paradox of Schrodinger's cat by experimentally con- 
firming that the cat was already dead before somebody 
opened the box to look. Likewise, we can conclude that 
the non-local collapse of entangled states merely corre- 
sponds to a Bayesian probability update in the remote 
system, providing a classical analogy that can explain 
non-classical properties of quantum mechanics such as 
the impossibility of non-local signaling and the transfer 
of quantum information by classical channels in quantum 
teleportation [HE!- 

VI. EXPERIMENTALLY OBSERVABLE 
NEGATIVE PROBABILITIES 

As explained above, weak measurement tomography 
identifies the quantum statistics of systems with well- 
defined initial and final properties. This means that the 
available information about each system can be more pre- 
cise than the uncertainty limit allows. Such super-certain 
information finds its quantum statistical expression in 
the non-positive density matrices Rif that summarize 
the results of all possible weak measurements between 
i and /. The negative eigenvalues of Rif represent weak 
values that exceed the eigenvalue bounds of positive den- 
sity matrices, providing a consistent framework for the 
resolution of quantum paradoxes by weak measurements 

Quantitative descriptions of quantum paradoxes such 
as Bell's inequalities [25[, Leggett-Garg inequalities [HH, 
6] or contextuality inequalities [26j are usually formulated 
in terms of precise measurements performed on different 
representatives of the same state. To connect this con- 
ventional formulation with weak measurements, it is nec- 
essary to show that the joint probabilities determined in 
weak measurements provide a unique and measurement 
independent definition of joint probabilities in quantum 
systems. In particular, the joint probabilities should not 
depend on the order in which the result is obtained [27| . 
For a pair of strong measurements {II/} and the 
joint probability from a weak measurement of g followed 
by a strong measurement of / should therefore be equal 
to the joint probability of a weak measurement of / fol- 
lowed by a strong measurement of g, 

P(f,9\i) = p(f\i)Tr{9 g Rif} 

= p(g\t)TT{fl f R lg }. (14) 

Weak measurement tomography can confirm that these 
two results are indeed equal. This means that the tran- 
sient density matrices Ri f and Ri g describe the same sta- 
tistical correlations between the measurement of / and 
the measurement of g. We can therefore conclude that 



weak measurement tomography provides a consistent def- 
inition of joint probabilities for measurements that can- 
not be performed jointly. In terms of the POVMs, this 
joint probability reads 

p(f,g\i) = Tt{~(* 9 n/ + nA) (15) 

It should be emphasized that each of these joint proba- 
bilities can be measured directly in an appropriate weak 
measurement. Eq. (|15p is therefore the only definition of 
joint probability with empirical validity. Interestingly, 
this means that quasi-probabilities such as the Wigner 
function should not be interpreted in terms of joint prob- 
abilities, since their mathematical construction does not 
refer to a valid pair of measurements 28] . On the other 
hand, the negative probabilities predicted by eq. (|15p are 
a natural consequence of quantum statistics, required by 
the consistency of weak and strong results. Thus, the va- 
lidity of the results is not a matter of interpretation, even 
though the precise meaning of negative joint probabilities 
for outcomes that never occur jointly may be difficult to 
understand. 



VII. IMPLICATIONS OF UNCERTAINTY FOR 
INDIVIDUAL QUANTUM SYSTEMS 

Obviously, negative probabilities cannot be interpreted 
as relative frequencies of actual measurements. Never- 
theless, they can be derived from the relative frequencies 
of weak measurements that consistently give the same 
results as the corresponding strong measurements. To 
reconcile the strangeness of negative probabilities with 
the empirical sense of reality justified by reproducible 
measurement results, it is important to remember that 
the validity of an anticipatory decomposition depends 
on the performance of the actual measurement. If an 
alternative measurement is performed instead, the sub- 
ensembles need to be "re-shuffled" before the correct de- 
composition is applied. For practical purposes, contra- 
dictions are avoided because the reality of an individual 
system is determined by only one of the possible measure- 
ments. The resolution of quantum paradoxes by negative 
probabilities is therefore based on the difference between 
the total statistical ensemble and its individual repre- 
sentatives. The reality of the representative is given by 
precise values of i and / (or i and g), while the ensemble 
properties are described by the density matrix. In the 
case of the sub-ensemble defined by i and /, this density 
matrix is given by Ri f , which may have negative eigenval- 
ues since it can only be observed in weak measurements. 
Eqs. (Til)]) show that the joint negative probabilities pre- 
dicted by weak measurement tomography of i and / are 
consistent with the joint negative probabilities obtained 
from i and g. The result can therefore be summarized 
in terms of a uniquely defined joint probability (| 1 5p that 
expresses the relation between strong measurements of 
/ and of g directly in terms of their POVM operators. 
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However, the negative eigenvalues of this operator clearly 
indicate that individual realities must be restricted to 
the actual measurement outcomes caused by the respec- 
tive system. Weak measurement tomography thus de- 
cides quantum paradoxes in favor of locality, causality, 
and non-contextuality, but against the notion of a non- 
empirical realism that attempts to provide a description 
of individual quantum objects beyond the uncertainty 
limited reality of its individual effects. 

VIII. CONCLUSIONS 

In conclusion, weak measurement tomography reveals 
a striking consistency of the quantum measurement for- 
malism with classical statistics by defining an unambigu- 
ous partition of the total ensemble described by fa into 
well-defined sub-ensembles Rif representing the future 
measurement outcomes /. The reconstruction of quan- 
tum states by weak measurements thus provides empir- 
ical evidence that the selection of a measurement out- 
come / does not eliminate the coherences between / and 
other outcomes. In particular, weak measurements can 
show that particles moving only through slit 1 of a double 
slit experiment carry the complete interference pattern of 
the initial state with them until the physical interaction 
of the final measurement randomizes the phase relation. 



Empirical evidence thus favors a statistical interpretation 
of quantum mechanics that assigns reality to individual 
measurement outcomes even before the measurement is 
performed. Quantum paradoxes can then be explained 
in terms of the negative probabilities described by the 
non-positive density matrices Rif. As was shown above, 
these transient density matrices uniquely define the joint 
probabilities between the measurement results / and the 
possible outcomes of other measurements g. The nega- 
tive values of such joint probabilities demonstrate that 
no consistent simultaneous assignment of both g and / is 
possible. Weak measurements can thus provide experi- 
mental proof that quantum paradoxes arise from the fal- 
lacious imposition of non-empirical realities beyond the 
uncertainty limits restricting the observable effects of in- 
dividual systems. 
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